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The magnetical field tuned superconductor-insulator transition in amorphous thin films, e.g., Ta 
and InO, exhibits a range of yet unexplained curious phenomena, such as a putative low- resistance 
metallic phase intervening the superconducting and the insulating phase, and a huge peak in the 
magnetoresistance at large magnetic field. Qualitatively, the phenomena can be explained equally 
well within several significantly different pictures, particularly the condensation of quantum vortex 
liquid, and the percolation of superconducting islands embedded in normal region. Recently, we 
proposed and analyzed a new measurement in Ref. 1 that should be able to decisively point to the 
correct picture: a drag resistance measurement in an amorphous thin-film bilayer setup. Neglecting 
interlayer tunneling, we found that the drag resistance within the vortex paradigm has opposite sign 
and is orders of magnitude larger than that in competing paradigms. For example, two identical 
films as in Ref. 2 with 25nm layer separation at 0.07K would produce a drag resistance ~ lO'^fi 
according the vortex theory, but only ~ 10^^^ SI for the percolation theory. We provide details of 
our theoretical analysis of the drag resistance within both paradigms, and report some new results 
as well. 



I. INTRODUCTION 

Amorphous thin film superconductors exhibit a vari- 
ety of fascinating quantum phenomena, due to the im- 
portance of fluctuation and disorder in two dimensions. 
Early theoretical"^ and experimental^"^^ work focus on 
the quantum superconductor-insulator-transition (SIT) 
in these materials. As one increases the perpendicular 
magnetic field or decreases the film thickness, the film 
changes from superconducting to insulating. An appeal- 
ing theoretical picture of the SIT is that the amplitude 
of the superconducting order parameter remains finite 
across the transition, and the transition is driven by 
phase fluctuations, which can be viewed as the conden- 
sation of vortices. Therefore the insulator is described 
as a vortex superfiuid, and the transition point is nearly 
self dual: it could be described either as the conden- 
sation of Cooper pairs, or of vortices. This Cooper- 
pair - vortex duality also suggests that the critical resis- 
tance at the transition should be _/?□ = h/Ae^ = 6.5kil, 
which is consistent with observations on strongly disor- 
dered samples"'^^. A variety of other experiments shows a 
transition with a critical resistance of the same order as 
Rq = h/Ae^. 

In recent years, experiments on these amorphous thin 
films have revealed more surprising results, mainly in 
transitions tuned by normal magnetic field. One of 
these raises the possibility that a metallic phase inter- 
venes between the superconducting and the insulating 
phases^'^°"^^. Near the "SIT critical point", as temper- 
ature is lowered below ^ lOOmK, the resistance curve 
starts to level off, indicating the existence of a novel 
metallic phase, with a distinct nonlinear I — V charac- 
teristics at least in Ta films that are interpreted as a 
consequence of vortex dynamics^'*. Another interesting 
experimental finding is the nonmonotonic behavior of the 



magnetoresistance^'^'^''^^'^®. As one increases magnetic 
field further from the " SIT point" , the resistance climbs 
up quickly to very large value in InO and TiN films, be- 
fore it plummeting back to the normal state resistance, 
as shown in FIG. 1. In Ta and MoGe films, as well as 
some InO films, the resistance peak is not as large, but 
is still apparent^^"^^ . 




FIG. 1: A typical magnetoresistance curve of amorphous thin 
film superconductors. As the magnetic field B increases, the 
superconducting phase is destroyed, and a possible metallic 
phase emerges. After which the system enters an insulating 
phase, where the magnetoresistance reaches its peak. The 
resistance drops down and approaches normal state value as 
B is further increased. 

Two competing paradigms may account for these phe- 
nomena. On the one hand, within the quantum vortex 



pictures 
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the insulating phase at the peak of the 



magnetoresistance implies the condensation of quantum 
vortices, and the high field negative magnetoresistance 
indicates the gradual depairing of Cooper pairs and the 
appearance of a finite electronic density of states at the 
Fermi level. The intervening metallic phase is described 
as a delocalzed but yet uncondensed diffusive vortex liq- 
uid as described in Ref. 31. In this picture disorder 



and charging effects are most important on length scales 
smaller or of order ^ (the superconducting coherence 
length, typically of order 10n?Ti). On the other hand, 
the percolation paradigm^^^^^ describes the amorphous 
film as a mixture of superconductor and normal or in- 
sulating puddles, with disorder playing a role at scales 
larger than ^. Particularly germane is the picture in Ref. 
35 which phenomenologically captures both a metallic 
phase as well as the strongly insulating phase by assum- 
ing superconducting islands exhibit a Coulomb blockade 
for electrons. This way the peak in the magnetoresis- 
tance arises from electron transport though the perco- 
lating normal regions consisting of narrow conduction 
channels. Yet a third theory tries to account for the low 
field superconductor-metal transition using a phase glass 
modeP^''^^ (see, however, Ref. 39 which argues against 
these results), but does not address the full magnetore- 
sistance curve. Qualitatively, both paradigms above are 
consistent with magnetoresistance observations, and re- 
cent tilted field'"^, AC conductance'*'^, Nernst effect"*^, 
and Scanning Tunneling Spectroscopic^^ measurements 
cannot distinguish between them. Particularly intrigu- 
ing is the origin of the metallic phase - is it vortex driven 
or does it occur due to electronic conduction channels 
dominating transport through the film? 
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FIG. 2: Our proposed bilayer setup for the drag resistance 
measurement. A current bias 7i is applied in one layer, and a 
voltage V2 is measured in the other layer. The drag resistance 
Rd is defined as Rd ~ V2//1. 



Given the similarity in the predictions of the dis- 
tinct vortex-condensation and percolation paradigms, an 
experiment that distinguishes between them would be 
highly desirable. We propose that a thin film "Giaever 
transformer"^'* experiment (FIG. 2) can qualitatively dis- 
tinguish between these two paradigms. The original de- 
sign of a Giaever transformer consists of two type-II su- 
perconductors separated by an insulating layer in perpen- 
dicular magnetic fields. A current in one layer moves the 
vortex lattice in the entire junction, yielding the same 
DC voltage in both layers. Determining the drag re- 
sistance Rd = V2//1 in a similar bilayer structure of 
two amorphous superconducting thin films should quali- 
tatively distinguish between the two paradigms (see also 
Refs. 45,46): within the vortex paradigm, vortices in 
one layer drag the vortices in the other, but within the 
percolation picture, the drag resistance is solely due to 



interlayer " Coulomb drag" , as studied in semiconductor 
heterostructures^^ . 

The first qualitative difference between vortex drag 
and Coulomb drag is the sign of the drag voltage V2. 
Denoting the voltage drop in the driving layer as Vi, it 
is easy to see that Vi and V2 have the same sign if they 
are produced by vortex motion, because vortices in the 
two layers move in the same direction transverse to the 
current bias /i. (We note in passing that if the second 
layer is in a closed circuit, the vortex drag would induce a 
current in the opposite direction in the secondary layer, 
since no outside voltage source balances the EMF pro- 
duced by the vortex motion.) On the other hand, Vi and 
V2 would have opposite signs if they are due to electron 
Coulomb drag, because V2 has to balance the drag force 
to ensure the open circuit condition in the second layer. 
In other words, Coulomb drag would try to produce cur- 
rent in the same direction in the primary and secondary 
layer. 

More importantly, we have found that in the vortex 
scenario, the drag resistance is expected to be several or- 
ders of magnitude larger than that in other models. Par- 
tially this is expected because in these films, the sheet 
carrier density ^ 10"'^^cm~^ is much larger than the vor- 
tex density ~ -B/$o ^ 10*°cm^^, and the drag effect 
is typically smaller for larger densities. For example, 
two identical films as in FIG. 2(b) of Ref. 2 with 25nm 
center-to-center layer separation at 0.07K would produce 
a drag resistance ^ 10~'*r2 according the vortex theory 
(see FIG. 3), but only ^ 10~*^ri for the percolation the- 
ory (see FIG. 4). But as we shall show below, the large 
vortex drag effect is also a consequence of the extremely 
high magneto-resistance slope, which has different impli- 
cations for the vortex condensation and percolation pic- 
tures. The strength of the thin-film Giaever tranformer 
experiment would therefore be in the transition region 
where the metallic phase transforms into the insulating 
phase, and the magneto-resistance is at a maximum. 

We believe that these qualitative differences between 
the drags in the two paradigms are quite general for 
each paradigm, and does not depend the various mi- 
croscopic assumptions made in various flavors of these 
phenomenological pictures. We will support these claims 
by analyzing the drag resistance between two identical 
thin films within a representative theoretical framework 
in the vortex"^* and percolation paradigms^*. We will re- 
strict ourselves to the standard drag measuring geometry 
assuming zero tunneling between the layers. We expect 
that allowing small tunneling will stregthen the effect; we 
will pursue this possiblity in future work. 

This paper is organized as follows. In Sec. II, we ex- 
tend the quantum vortex formalism to bilayers, and then 
we calculate the drag resistance in the insulating and the 
metallic regime, respectively. The effect of unpaired elec- 
trons on the drag resistance is also studied. In Sec. Ill, 
we review the percolation theory of Ref. 34, and then ex- 
tend this theory to bilayers as well, in order to calculate 
the drag resistance. In Sec. IV, we briefiy discuss the 



drag resistance behavior within the phase glass model of 
Refs. 37,38. Finally, we summarize and discuss our re- 
sults in Sec.V. Some details are provided in appendices. 



II. DRAG RESISTANCE IN THE QUANTUM 
VORTEX PARADIGM 



A. The vortex description of double-layer 
amorphous films 

Within the quantum vortex paradigm, the insulating 
phase has been explained as a superfluid of vortices by 
the "dirty boson" model of Ref. 4, while the metallic 
phase is expected to be an uncondensed vortex liquid (see 
also Ref. 29). This picture has been pursued by Ref. 31 
which argues that vortices form a Fermi liquid for a range 
of magnetic field, thereby explaining the metallic phase. 
At larger fields, where the insulating phase breaks down, 
it is claimed that gapless bogolubov quasi particles nick- 
named spinous, i.e., unpaired fermions with finite density 
of states at the Fermi energy, become mobile, impede vor- 
tex motions, destroy the insulating phase, and suppress 
the resistance down to normal metallic values. 

We will concentrate on the case where no interlayer 
Josephson coupling exists, and the vortex drag comes 
from the magnetic coupling between vortices in different 
layers which tends to align themselves vertically to min- 
imize the magnetic energy. To calculate the drag resis- 
tance in a bilayer setup, it is crucial to derive the vortex 
interaction potential due to the current-current magnetic 
coupling between the layers, which is captured by the B^ 
term in the Maxwell action. We achieve this by both a 
field theory formalism and a classical calculation. The 
classical calculation is relegated to Appendix C. 

Let us next derive the vortex action. Treating the su- 
perconducting film as a Cooper pair liquid, we have the 
following partition function 



Z = Vp{Dp2Ve{De2VAe 
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where a is the (center-to-center) layer-separation, p„ 
and 6n are the 2d density and phase fiuctuation of the 



ri— th layer Cooper pair field, respectively, A and A^xt 
are the fluctuating and external part of the electromag- 
netic field, respectively. The intralayer Coulomb interac- 
tion Vi{r) = (2e)^/r (whose 2d Fourier transform would 
be 27r(2e)^/g), and the interlayer Coulomb interaction 
Ve{r) = {2e)^ / y/r^ + a^ (whose 2d Fourier transform is 
27r(2e)^/(7e^''"). ps is the superfiuid phase stiffness of 
each layer, which can be determined approximately from 
the Kosterlitz-Thouless temperature Tkt' 



Tkt ^ ^Ps 



(2) 



Next, we follow a procedure of vortex-boson duality 
transformation taking into account the B^ term (which 
will be the origin of the interlayer vortex interaction), 
and obtain the following dual action for the vortex field 
i/i^jn of the n-th layer and two U(l) gauge fields a^ and 
/3^ (see Appendix B for details): 
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(3) 



where Sp^n — Pvn ^ B/^q, $o is the flux quantum, p-^n = 
tplj^tpvn, <t>n is the phase of the vortex field ipyn, and my is 
the vortex mass. Since there is still controversy over the 
theoretical value of m„, we chose to determine the vortex 
mass from experiments. As discussed in Appendix A, for 
the InO film of Ref. 2, we obtain rriy « 19?ne where We 
is the bare electron mass. 

a^ and I3i, are gauge fields which mediate the symmet- 
ric and antisymmetric part of the vortex-vortex interac- 
tion. They are related to the Cooper pair currents j„^ in 
the n— the layer by 
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(4) 



For i^ = 1,2, the dual charges and the dual "light 
speeds" are 
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(6) 



where qc is the inverse of the 2d Pearl screening length^^, 
which can be estimated from the value of Tkj-'- 
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For example, the film in Ref. 2 has Txt around 0.5K. 
This corresponds to qc ~ (4cm)~^, and it is much smaller 
than the inverse of typical sample size 1/L ^Imm"'^. 

In (3), we have chosen the transverse gauge for the 
gauge fields a^ and /S^ and integrated out ao and /3o to 
obtain the vortex interaction potentials. The intralayer 
vortex interaction potential 



U^{q) 



nic 
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and the interlayer vortex interaction potential 



Ue{q) 



-""c/. 



(8) 
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When r < 1/qc, Ui{r) gives the familiar log interaction; 
for r > 1/qc, i-C, beyond the Pearl screening length, 
Ui{r) is still logarithmic but with half of the magnitude^'', 
in contrast to the l/r behavior of the single layer case 
(which is Eq. (8) with a — > oo). The interlayer inter- 
action Ue is purely due to the magnetic coupling, i.e., 
vortices in different layers tend to align to minimize the 
energy cost in the B^ term. As expected, the interaction 
between two vortices with the same vorticity in different 
layers is attractive, although its strength is suppressed 
with increasing distance a and decreasing qc. Ui and Ue 
can also be derived classically by solving London equa- 
tions and Maxwell's equations, which we will show in 
Appendix C. In addition, the form of Ue is equivalent to 
those derived in Ref. 50,51. 

Following Ref. 29, one can examine the strength of the 
interaction between vortices and transverse gauge field 
modes by looking at the dimensionless coupling constant 
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for the entire range < q < 1/^, ^ ^ lOnm being the 
coherence length. Thus, the transverse gauge field exci- 
tations can be neglected. For a comparison, the dimen- 
sionless parameter for the strength of the longitudinal 
interactions Ui and Ue is 



,m,; 
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(11) 

With these simplification, we now rewrite the action 
for the bilayer system as 
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As the magnetic fields increases, Ol gets suppressed, 
and therefore the vortex system goes from a interaction- 
dominated localized phase (Cooper-pair superfluid phase, 
i.e., superconducting) to a kinetic-energy-dominated su- 
perfluid phase (Cooper-pair insulating phase), possibly 
through a metallic phase. Finally, when the applied 
magnetic field is large enough that unpaired electrons 
("spinous" in Ref. 31) are delocalized, they impede vor- 
tex motion through their statistical interaction with vor- 
tices and therefore suppress the resistance down to values 
consistent with a normal state in the absence of pairing 
(see Ref. 31). 



B. Drag resistance in the vortex metal regime 

As explained in the introduction, essentially all films 
undergoing a magnetic field driven SIT also exhibit the 
saturation of their resistance at the transition. Within 
the vortex picture, the intervening metallic phase is in- 
terpreted as a liquid of uncondensed vortices"^^, and the 
vortices are diffusive, and have dissipative dynamics. At 
intermediate fields and low temperatures, where the in- 
termediate metallic phase appears, the vortices are delo- 
calized but uncondensed. In this phase one can derive the 
following form of the the drag conductance an (which for 
the vortices is the equivalent through duality to the drag 
resistance of charges) using either the Boltzman equation 
or diagrammatic techniques, irrespective of the effective 
statistics of vortices'*^'^^"^^: 

h^ daidn2 r 3 a Ha I^^P ImXi IniX2 
"- = 8^^a;^X '^'^Vo '" sinh^(l^) ' 

(13) 
where Ui, rii, and Xi ^'^^ the conductance, density, and 
the density response function of the vortices in the i— th 
layer. In addition. 



U = 



Ue 
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is the screened interlayer interaction, Ue is the bare inter- 
layer interaction, and Ui is the intralayer interaction, and 
T is the temperature, da^/dn^ appears since Rb is re- 
lated to the single layer rectification function, F, defined 
^s jy = Tcjy^ , with (j> being the vortex potential field. F is 
generally proportional to dav/duv (see Ref. 56). Com- 
bining the vortex density expression Ui = B/^q and the 
relation between physical resistance and the vortex con- 
ductance R = (j-)^crt, with (13), one obtains the drag 
resistance 



Rd 



e^% dRi dR2 

Stt^T dB dB 



q^dq / duj\U\ 
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sinh^ (If) 



(15) 

Remarkably, the drag resistance is proportional to 
dRi.2/dB, and thus Rjj peaks when the MR attains its 
biggest slope. This is one of the most important results of 
our analysis. Intuitively, the dependence of the drag on 



day /dnv — dRi,2/dB arises since the drag effect is the 
resuh of the nonuniformity of the relevant particle den- 
sity; how this nonuniformity affects the voltage drop in 
the medium both in the primary and secondary layers is 
exactly the origin of the square of the magneto- resistance 
slope. 

The only model-dependent input is the density re- 
sponse function Xi.2- We have computed the drag re- 
sistance using two different choices of Xi,2- In the re- 
mainder of this section, we follow the vortex Fermi liquid 
description for the metallic phase of Ref. 31 and use the 
fermionic response function for X1.2; in Appendix D, we 
treat the metallic phase as a classical hard-disk liquid 
of vortices^^'^^ and use its response function accordingly 
for Xi,2- It turns out that the drag resistance results are 
remarkbly close for these two approaches, hence showing 
the robustness of our results. 

If wc treat vortices as fermions in this phase'^^, we use 
the Hubbard approximation form for xi.2 considering the 
short-range repulsion between vortices and also the low 
density of this vortex Fermi liquid^^'^'': 



xiq,^) = 



xo(g,w) 



l-U^{q)xo{q,Lo)G{q)^' 



(16) 



where G{q) — q^ /{q^ + kp), and kp of the vortex Fermi 
liquid can be easily calculated from the vortex density: 
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One can define the mean free path / and the transport col- 
lision time T for vortex Fermi liquid. Their value can be 
estimated by combining the expression for vortex conduc- 
tivity ay — nyT/niy and the relation between the physical 



resistance and the vortex conductance R = (^) f 




(18) 



When ql > I or lot > 1 we approximate xo by the non- 
interacting ballistic fermion result^^: 
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where 
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and 
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(21) 



For ql < 1 and lot < 1, we use the diffusive Fermi liquid 
result: 



Xa = v 



Dq^ 



Dq^ — iu) 



(22) 



Plugging (16) into (15), one can numerically compute the 
drag resistance. The result is given in Sec. HE. 

Note that this result does not crucially depend on 
choice of fermionic density response function above. As 
stated earlier, as long as vortices form an uncondensed 
liquid, (15) remains valid. We have also computed Rd 
by modeling the metallic phase as a classical hard-disk 
liquid of vortices'"'^'^^, and putting the corresponding den- 
sity response function into (15). The resulting magnitude 
and the behavior of Rjj are extremely close to the re- 
sults we obtained above within the vortex Fermi liquid 
frameworks (see Appendix D). This demonstrates the 
universality of our results. 



C. Drag resistance in the insulating (vortex 
superfluid) regime 

According to the vortex theory, the insulating phase is 
a superfluid of bosonic vortices. In this regime, the vortex 
dynamics is presumably nondissipative. A mechanism of 
nondisspative supercurrent drag between bilayer bosonic 
superfluid systems has been studied by Ref. 62-64. Here, 
we apply this approach to the superfluid of vortices in 
the insulating regime. In the absence of current bias, we 
have the following action from (12) deep in the insulating 
phase: 



s = J2{-'^p 



2,m„ 



- iSp2(j>20J + 7r^i-q^<p2) 

+ \u^(5plf + \ud5p2f + Ue5pi5p2\ . (23) 

Switching to the canonical quantization formalism and 
using mean field approximation for the quartic interac- 
tion term^^, the above action (23) corresponds to the 
following Hamiltonian for bilayer interacting bosons: 

^ = E T.[4;r''^s{q)as{q) + ^[UM) + SUM] 



X H{q)al{-q) + as{-q)a,{q)]} 



where 



a±{q) = -^[tpviiq) ± ■4'v2{q)] 



(24) 



(25) 



'011 1 and '01)2 are the bosonic vortex field operators for the 
first and second layer, respectively. (24) can be diagonal- 
ized using Bogoliubov transformations: 

a±{q) = M±(<?1&±(<?1 + v±{q)bi{-q), (26) 



where in the long wavelength limit 
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A vortex current bias vl in layer 1 (the driving layer) is 
represented by a perturbation term Hi in our Hamilto- 
nian: 



Hi =^ myji 



■Vl. 



(28) 



The drag current in the second layer can be calculated 
using standard perturbation theory. The new ground 
state to the first order in vi is given by 
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where |0) is the vacuum state of 6lj_, and \n) represents 
all possible states obtained by acting 6]j_ on |0). Thus, at 
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It is straightforward to check that the only excited 
states \n) that contribute to the sum are of the form 
b_^{q)b_{—q)\0). One thus obtains 
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(31) 



Now, plugging (27) into (31), to the second order in in- 
terlayer interaction Ue we have 



in) = vl 
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D. The effect of mobile spinons 

The discussions in previous sections apply to the case 
where no mobile unpaired electrons, i.e. spinons in Ref. 
31, exist in the system. However when the magnetic field 
is strong enough to pull apart Cooper pairs and delocalize 
spinons, as is signaled by the downturn of the niagnetore- 
sistance, the drag resistance is modified by the spinons. 
In this subsection, we analyze how mobile spinons affect 
our drag resistance results above. 

We follow the semiclassical Drude formalism as in Ref. 
31 which takes into account the statistical interaction 
between Cooper pairs, vortices, and spinons. Vortices 
and spinons see each other as 7r-flux source, while electric 
current exerts Magnus force on vortices. Denoting the 
electric current, vortex current, and the spinon current 
in the n— th layer as J„, jv,m js,n, we have the following 
equations for the first (driving) layer (see Ref. 31): 

jyi ^ ayzx {jsi - Ji), 

jsl ^ (TsZX jyl. 

Similarly, denoting the vortex drag conductance without 
spinons as ctd, we incorporate the drag effect in the fol- 
lowing way in the equations of the second (passive) layer: 

Jd2 = .?t>l + CTyZ X Js2, 

C7y 

js2 ^ (TsZ X jy2- 

This set of equations is a consequence of the absence of 
electric current but the presence of vortex drag effect in 
the second layer. We can solve these two sets of equa- 
tions, and obtain the effective vortex drag conductance: 
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(33) 



Since the physical resistance R = {h/{2e))'^ay, we have 



R 



<^ff 



Rt 



(1 + Ry/R, 



(34) 



where Rd is the drag resistance if spinons are localized, 
Ry = {h/2e)'^(jy is the vortex contribution to the resis- 
tance, and Rs = oj^ is the spinon contribution to the 
resistance. Thus, we see that when Rs ^ Ry, the drag 
resistance is quickly suppressed to unmeasurably small 
as spinon mobility increases. 



Divding this result by (ji) = riyVi and recalling that the 
resistance is proportional to the vortex current, one is 
ready to obtain the drag resistance. 
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(32) 



When spinons are mobile, they will suppress the drag 
resistance, as we will show in section II D. 



E. Results of the drag resistance in the vortex 
theory 

Collecting the above results and the value of the vor- 
tex mass ruy discussed in Appendix II A, tuning the value 
of the vortex (spinon) contributions to the resistance Ry 
(Rs) so that R = RyRg/{Ry + Rg) (see Ref. 31) resem- 
bles the resistance observed in the experiment of Ref. 2, 
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FIG. 3: Drag resistance Rd (in Ohms) between two identical 
films as in FIG. 2b of Ref. 2 vs. magnetic field B, accord- 
ing to the vortex picture^^ (log scale); . The drag resistance 
has been smoothened to avoid discontinuity at the boundary 
between the metallic and the insulating phase. Center-to- 
center layer separation a — 25nm, temperature T = 0.07K 
and 0.35K. Insets: single layer magnetoresistance (magne- 
toresistance, log scale) reproduced according to the quantum 
vortex theory.. The parameters are tuned to make the mag- 
netoresistance resemble the experimental data in FIG. 2b of 
Ref. 2. Rd has a peak at the steepest point (~ 8T) of the 
magnetoresistance, which is due to the fact that Rd is pro- 
portional to the square of the slope of the magnetoresistance 
in the small magnetic field side of the peak. Also, Rd is 
larger at lower temperature, because the magnetoresistance 
curve is then much steeper. Carrying out the experiments at 
even lower temperatures may further enhance the vortex drag 
effect. 



and setting temperature to be 0.07K and 0.35K, we have 
calculated the drag resistance between two identical films 
with single layer resistance given by the inset of FIG. 3, 
and with centcr-to-center layer separation 25nm. We as- 
sume that vortices form a Fermi liquid (thus (15) is appli- 
cable; however see also Appendix D) when B < 9T, and 
a bosonic superfiuid (thus (32) is used) when B > 9T. 
We smoothcn the drag resistance curve by convoluting 
it with a Gaussian function to avoid discontinuity across 
the phase boundary between the metallic phase and the 
insulating phase. 

The results of vortex drag are summarized in FIG. 3. 
One can see that The drag resistance has a peak at the 
steepest point (~ 8T) of the magnetoresistance. This 
is due to the fact that in the vortex metal regime, the 
drag resistance is proportional to the square of the slope 
of the magnetoresistance. Also, the drag resistance is 
larger at lower temperature. This is because the mag- 
netoresistance curve is much steeper as one approaches 
zero temperature (see (15)). For the film of Ref. 2, the 
sheet drag resistance is about 10~^ mil at its maximum, 
which is measurable despite challenging. We suggest to 
carry out experiments to even lower temperature, which 
should leads to a larger drag resistance. Using a Hall-bar 



shape sample would also amplify the result. 



III. DRAG RESISTANCE IN THE 
PERCOLATION PICTURE 

A. Review of the percolation picture of the 
magnetoreistance 

Within the percolation picture of Ref. 34, it is ar- 
gued that the non-monotonic magnetoresistance arises 
from the film breaking down to superconducting and nor- 
mal regions (described as localized electron glass) ^'*. As 
the magnetic field increases, the superconducting region 
shrinks, and a percolation transition occurs. Once the 
normal regions percolate, electrons must try to enter a 
superconducting island in pairs, and therefore encounter 
a large Coulomb blockade absent in normal puddles. The 
magnetoresistance peak thus reflect the competition be- 
tween electron transport though narrow normal regions, 
and the tunneling through superconducting islands. 

This picture is captured using a resistor network de- 
scription. Each site of the network has a probabil- 
ity p to be normal, and 1 — p to be superconducting; 
each link is assigned a resistance from the three values 
Rnn, Rss, Rsn, that reflect whether the sites the link 
connects are normal (N), or superconducting (S). An in- 
crease of the magnetic fleld is assumed to only cause p 
to increase. Since the normal region is described as dis- 
ordered electron glass, Rnn, the resistance between two 
normal sites, is assumed to be of the form of hopping 
conduction: 
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where ^loc is the localization length, and e^ is the energy 
of the I— th site measured from the chemical potential 
(taken from a uniform distribution [—W/2,W/2]), and 
for simplicity we allow only nearest neighbor hopping. 
The resistance between two superconducting sites, Rss, 
is taken to be very small, but still nonzero, and vanishes 
as T ^ T" —> 0. Most importantly, the resistance be- 
tween one normal site and a neighboring superconducting 
site, RsN, is assumed activated: 



RsN ^ RsNo exp 



knT 



(36) 



to model the charging energy electrons need to pay to 
enter a superconducting island. 

We have reproduced the work of Ref. 34 where the 
parameters of this model are chosen to reproduce the 
magneto-resistance curves and temperature dependence 
observed in the strong- insulator InO sample^. The total 
resistance vs. the probability of normal metal (assumed 
to increase with increasing magnetic fleld) is shown in 
the inset of FIG. 4. Indeed, the peak of the magnetore- 
sistance can be explained by this theory. However, as we 



demonstrate now, this theory predicts a very different 
behavior for the drag resistance. 



Calculation of drag resistance within the 
percolation picture 
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FIG. 4: Drag resistance Rd (in Ohms) between two identical 
films as in FIG. 2b of Ref. 2 vs. normal metal percent- 
age p (corresponding to normal magnetic field) , according to 
the percolation picture'^'*. Center-to-center layer separation 
a = 25nm, temperature T — 0.07K and 0.35K. Insets: single 
layer magnetoresistance (magnetoresistance, log scale) repro- 
duced according to the percolation theory. The parameters 
are tuned to make the magnetoresistance resemble the experi- 
mental data in FIG. 2b of Ref. 2. The sign of the voltage drop 
of the passive layer is opposite to that of the driving layer, 
and the maximum magnitude value of Rd is much smaller. 

To calculate Rjj, we first follow Ref. 34 and tune the 
parameters to make the single layer resistance resemble 
the experimental data in FIG. 2(b) of Ref. 2: £^ioc — 
0.1, W ^ 0.4K, Ec = 0.6K, Rsm - 10^^, and Rm - 
10~^r2. Next, we place one such network (active layer) 
on top of another one (passive layer) . Each link is treated 
as a subsystem, which might induce a drag voltage (an 
emf) £ = IRd in the link under it in the passive layer. 
When a link is between two normal (or superconducting) 
sites, it is treated as a disorder localized electron glass (or 
superconductor). In Appendix E, we find Rd between 
two localized electron glass separated by vacuum is: 
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Here, n k, b x lO^'^cm^'^ is the typical carrier density of 
InO^"^, d ~ 20nm is the film thickness, a = 25nm is the 
center-to-center layer separation, Ri^2 are the resistances 
of the two normal-normal(NN) links, xo — a/{2TTe'^vd^'^) 
where v is the density of states and ^ « Inm is the lo- 
calization length. The value of the localization length 



^ is estimated by following Ref. 34 to take ^ ~ 0.1 x 
plaquette size (reflecting the fact that it is a disordered 
insulator), and we estimate the plaquette size as the su- 
perconducting coherence length ^ lOnm. Although this 
estimation of localization length is crude, the drag resis- 
tance Rd has only logarithmic dependence on it in (37). 
Setting T = 0.07K, and R^ ^ R2 ^ lO^fl, we can esti- 
mate Rd ^ 10^12^. 

On the other hand, we will show in Appendix F that 
a genuine (i.e., without mobile vortices) superconductor 
has no drag effect at all in a resistor network, either 
when it is aligned with another superconductor link or 
a normal link. Thus, drag effects associated with a su- 
perconducting link can only come from vortices. How- 
ever, The small resistance for the superconducting is- 
lands in this theory implies that vortices in the super- 
conducting islands, if any, have very low mobility. If 
two superconducting links are vertically aligned, we can 
estimate the drag resistance due to mobile vortices us- 
ing our vortex drag result (15): roughly Rd oc i?^, 
for R ^ lO^ri we obtained Rd ^ 10~^r2, therefore for 
i? ^ 117 we have Rd ^ lO^^'^fl, which is negligible com- 
pared to the Coulomb drag resistance between two NN 
links ^ lO^^^il. Finally, Ref. 65 has shown that a cur- 
rent off the plane where vortices reside does not exert any 
force on vortices. By Newton's third law or equivalently 
the Kubo formula for the drag conductance, this also im- 
plies that moving vortices does not exert any DC emf in 
another layer. Therefore, there is no drag effect when 
a NN link is aligned with a SS link. Consequently, the 
Coulomb drag between two vertically aligned NN links 
(Eqn. (37)) dominates the drag effect. 

Thus, we solve the Kirchoff's equations for the two 
layers, and obtain the voltage drop and thereby the 
drag resistance. The results are shown in FIG. 4, with 
T = 0.07K and 0.35K, film-thickness 20nm, and the 
center-to-center interlayer distance 25nm. We observe 
that the sign of the voltage drop of the passive layer is 
opposite to that of the driving layer (not shown in the 
Figure), as expected and explained in the introduction, 
and the maximum magnitude of the drag resistance is 
around 10~^^il, indeed much smaller than that in the 
vortex paradigm. 



IV. DISCUSSION ON THE DRAG RESISTANCE 
IN THE PHASE GLASS THEORY 



A third theory, namely the phase glass theory^^''^®, fo- 
cuses on the nature of the metallic phase intervening the 
superconducting and insulating state. In this theory, the 
system is described as interacting bosons (Cooper pairs), 
but it is argued that the glassy phase is in fact a Bose 
metal, due to the coupling to the glassy landscape. 

Specifically, Ref. 37 has studied the quantum rotor 



model 
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where the Josephson couphng Jij obeys a Gaussian dis- 
tribution with nonzero mean. This model is appears 
to exhibit three phases: superconducting phase, phase 
glass phase, and a Mott insulator phase. Ref. 37 has 
employed replica trick to obtain the Landau theory of 
the the phase glass phase near the glass-superconductor- 
transition critical point, and has calculated the conduc- 
tance in this regime. It was found that in this regime the 
DC conductance is actually finite at zero temperature. 
For completeness, we note that Ref. 39 argued against 
these results and obtained infinite conductance instead. 
This analysis has recently been extended to include the 
external perpendicular magnetic field"^®, which is more 
relevant to the experiments on the magnetic field tuned 
transition. However, Ref.38 has only studied the regime 
of small magnetic field where one just enters the resistive 
glassy phase and left out issues such as the peak in the 
magnetoresistance. Therefore, we leave a complete anal- 
ysis to future work and simply observe that according 
to this theory, the resistive state is a glassy phase where 
phase variables 6i 's of the bosons are ordered locally. In 
other words, there are no mobile vortices moving around. 
Consequently, the current coupling as we considered in 
the vortex drag should is absent, and the Coulomb in- 
teraction should dominate the drag effect. Therefore, we 
expect that the sign of the drag voltage is opposite to 
the voltage drop of the driving layer, as we discussed in 
the introduction to be a general feature of the Coulomb 
drag, and the magnitude of the drag resistance should be 
small. This is in part because for a bosonic system, the 
phase space available for excitations is much smaller than 
fermionic systems due to the absence of a Fermi surface. 



V. SUMMARY AND DISCUSSION 

One of the most exciting possiblities is that the SIT 
in amorphous thin films realizes the vortex condensation 
scenario^'^'^^. The amorphous-films Giaver transformer 
experiment^, would be able to measure a distinct signa- 
ture of mobile vortices, which is a drag resistance oppo- 
site in its direction to that of coulomb drag. Therefore 
such a measurement would able to disclose whether the 
vortex paradigm is suitable for explaining the complex 
phase diagram of amorphous films in a normal manetic 
field, or whether the percolation paradigm is indeed more 
appropriate. We provide a detailed computation of the 
drag resistance according to the vortex theories of Ref. 
4,31 and the percolation theory of Ref. 34. The drag 
resistance implied by the phase glass model^^'^^ is also 
briefly discussed. We find that vortex picture predicts a 
drag resistance orders of magnitude stronger than non- 
vortex pictures. In addition, the drag resistance and the 



single layer resistance have the same sign according to 
the vortex picture, but the opposite sign for non-vortex 
pictures. Therefore, drag resistance measurement are in- 
deed able to distinguish different theoretical paradigms 
qualitatively. 

Wc considered specifically a bilayer device which will 
contain two identical films as in Ref. 2 with 25nm layer 
separation and at 0.07K. A calculation within the vor- 
tex paradigm yields a drag resistance Rd ^ lO^^il at 
its maximum value. This drag arises solely from the 
attractive interaction of the demagnetizing currents of 
vortices. The value we find is probably near the limit 
of measurability; we suggest, however, to carry out ex- 
periments at even lower temperature, in which case the 
single layer magnetoresistance is even steeper, and the 
drag resistance should be larger. Within the percola- 
tion picture of Ref. 34, the dominating drag effect is 
the drag between two vertically aligned normal regions 
in the different layers. For two identical films as in Ref. 
2 with 25nm layer separation at 0.07K, we find the drag 
resistance Rjj ^ 10~^^il at its maximum value, which is 
indeed orders of magnitude smaller than the drag resis- 
tance predicted by the vortex picture. Also, we find the 
sign of the drag resistance is the opposite of that of the 
single layer resistance, as expected. 

The answer we find should not depend crucially on 
the details of the microscopic picture which we use. If 
vortices are not responsible for the inhibitive resistance 
which the films display, then drag effects will appear 
primarlily due to Coulomb repulsion of single electrons. 
This drag effect will be low because of the relatively high 
electronic density in the films. On the other hands, if vor- 
tices are responsible for the large resistance in the inter- 
mediate magnetic fields leading to the insulating phase, 
then they will produce a drag opposite in its direction 
to the Coulomb drag. To carry out the vortex drag cal- 
culation in the metallic phase intervening between the 
superconducting and insulating phase we used the pic- 
ture of Ref. 31, which treats the vortices as fermionic 
diffusive particles. This picture is justified due to the 
strong long-ranged interactions within the vortex liquid, 
which render the question of statistics secondary, intu- 
itively, since vortices rarely encircle each other. Never- 
theless, to demonstrate the universality of our results, we 
also carried out the drag calculation in the metalic phase 
assuming that the vortices are hard core disks, and ob- 
tained essentially the same answer (c.f. App. D). 

Indeed our strongest results are obtained in the 
intermediate-field metallic phase. The controversy sur- 
rounding this phase requires some special attention. 
First, we note that all experiments of thin amorphous 
films exhibit a saturation of the resistance at temper- 
ature below about lOOmK at intermediate resistances. 
This is clearly seen in, e.g, the resistance vs. field traces 
which overlap at subsequent temperature sweeps as in 
Fig. 2b of Ref. 2. Second, there are reasons to believe 
that this saturation is not the result of failure to cool 
electrons. Resistances that are too low or too high con- 
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tinue to change as the temperature is lowered. But the 
two heating mechanisms most hkely are current heat- 
ing, with power ^ PR, and therefore affecting the high- 
est temperatures, and ambient RF heating, which would 
have a voltage-biased power ^ V'^/R, and therefore most 
effective in the lowest resistances. Neither mechanism 
explains resistance saturation at intermediate tempera- 
tures. Furthermore, experiments on Tantalum films show 
distinct signatures in the metallic regime which disap- 
pear in the insulating and superconducting regimes, and 
also distinguish it from the thermally-destroyed super- 
conducting phase^**. Third, even if the metallic behavior 
of the films is a finite temperature phenomena, within 
the vortex paradigm, the resistance still arises due to 
vortex motion. Therefore the drag calculated within this 
paradigm using a diffusive vortex model should still be 
adequate, and our results do not depend crucially on 
the existence of a zero-temperature intervening metallic 
state. 

The signatures we expect to find in the proposed mag- 
netic and Coulomb drag measurements are not large. In- 
corporating interlayer electron and Josephson tunneling 
will increase both the vortex-drag effect and the com- 
peting Coloumb drag effects. As we point out here, the 
drag signature of vortex motion, or single electrons or 
Cooper-pairs motion will have opposite signs. Quite pos- 
sibly, allowing interlayer tunneling will render both drag 
effects measurable. Indeed, such a setup will be a de- 
viation from standard drag measurements where charge 
transfer between layers is forbidden. Nevertheless, a care- 
ful choice of tunneling strength and sample geometry will 
make such experiments plausible and useful. We intend 
to analyze the vortex and Coloumb drag in the presence 
of interlayer tunneling in future work. 



where the action S is 
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Here, p and 9 are the density and phase fluctuation of 
the Cooper pair field, respectively, A is the fluctuat- 
ing electromagnetic field, and Apxt is the applied exter- 
nal electromangetic field, typically a perpendicular mag- 
netic field. V{r) = (2e)^/r (whose 2d Fourier transform 
would be 2ii{2ef' /k) is the Coulomb interaction between 
Cooper pairs, ps is the bare stiffness for phase fluctua- 
tions. The value of pa can be determined approximately 
by the zero-field Kostcrlitz-Thouless temperature Tkt'- 



Tkt = -^Ps- 



The 2d number current of Cooper pairs is 
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One can introduce the dynamical field j by Hubbard- 
Stratonavich transformation (or Villain transformation 
in the lattice version of this derivation) and transform Z 
to be 
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Appendix A: The determination of the vortex mass 

In this appendix, we demonstrate in detail the deriva- 
tion of the vortex-boson duality for a single layer and 
discuss the value of the vortex mass. Our starting point 
is the following partition function for Cooper pairs: 



Here, i is the imaginary number unit, q is the in-plane 
2d wave vector, while kz is the 3rd wave vector compo- 
nent perpendicular to the plane, and subscripts q mean 
Fourier transformed variables. Next we split the 6 field 
into a smooth part 6s and a vortex part 9y: — Og + Oy. 
Afterwards one can integrate out 9s to obtain the conti- 
nuity constraint: 
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Integrating by parts, and noting the definition of the vor- 
tex current density 
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we obtain 
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Furthermore, noting that A{q, z = 0) = J ^A{q,kz), 

one can integrate out A in its transverse gauge, and the 
action S now reads 
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where $0 — hc/{2e), and the "dual charge" of vortices is 

(A17) 



(A8) 



where qc is the inverse of the 2d Pearl screening length^^, 
and typically it is much smaller than 1/ L, where L is the 
sample size. 

The continuity constraint is solved by defining a new 
gauge field a^ — (oq, a) such that 
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where j^ — {c*p,j) and 9^ — (^Q^., V), and the value of 
constant 77 and the "speed of light" c* are to be deter- 
mined. Writing in components, (A9) is 
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where e and h are the dual " electric field" and " magnetic 
field" associated with a, respectively. To fix -q and c*, we 
require 

thus 
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Using (A9), we express the partition function Z as 

Z = fvaVaoVeye-'^, (A13) 
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In the above, we have assumed that the only external 
electromagnetic field is a perpendicular magnetic field 

Bext- 

The magnitude of the Magnus force, which now ap- 
pears as the electric force, can be easily verified: 

F = e* x\e\ = ^^ XTjj = hj, (A18) 

77 

as expected. 

Introducing a vortex field ip^ and making the action 
explicitly gauge- invariant, we write the action as 
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where ip^, — Pv — j^ , and we have introduced the vortex 
mass niy. Integrating out oq, one obtains 

S = y ^ <^ —Spyhiujcf) + -6pvUSpv 
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is the well-known Pearl interaction potcntial^^. 

In the insulating phase, i.e., the vortex condensed 
phase with vortex superfluid stiffness p^s, we have 
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Due to the Higgs mechanisni in this "symmetry broken 
phase" , the gap of the two modes in the vortex superfluid 
phase coincide to be 



Egap = yj2T:p^sel 
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for qc <C L^^ . Roughly speaking the two modes corre- 
spond to a density fluctuation of the vortices, or of the 
underlying Cooper-pairs Deep in the insulating phase, 
i.e., near the peak of the magnetoresistance, the vortex 
stiffness is simply 



Pvs 



(A24) 



where the vortex density n^ — -B/$o- Therefore, in this 
regime we have 



Egap = 2nh. ps 



(A25) 



Since the gauge field a^ is actually the fluctuation of 
Cooper pairs, we conjecture that its gap Egap can be 
identified with the activation gap observed in the exper- 
iments of Rcf. 2,23 near the insulating peak. Ref.2,23 
have also found that with increasing disorder strength, 
the ratio Egap/TuT is enhanced. This is natural from 
our expression (A25): dividing (A25) by (A3), we have 
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increasing disorder makes vortices more mobile and 
thereby suppresses the vortex mass mj'\ it also sup- 
presses the superfluid stiffness ps- Therefore, Egap /Tkt 
is larger for more disordered sample. 

Since there is still controversy over its theoretical value, 
we chose to use the experimental value of Egap as an 
input to deduce the vortex mass from (A25). Combining 
(A3), we can express the vortex mass m^ as a function 
of observable quantities: 



m^ 



~^ ■ 

gap 



(A27) 



Again, the vortex density n^ — S/$o- For the InO film 
of Ref. 2, Tkt ~ 0.5K, and Egap ~ 1.6K at B = 9T. 
Plugging these into (A27), we obtain my « 19me where 
TOe is the bare electron mass. For comparison, this value 
is not far from that of the so-called core mass of dirty 



superconductors' 
carrier density 
2,23). 
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m ~ {kFd)me 



49me if we use 
5 X 10^°cm~^ and d ~ 20nm (see Ref. 



Appendix B: The field theory derivation of the 
vortex interaction potentials in bilayers 

For identical bilayer superconducting thin films sep- 
arated by a (center-to-center) distance a, we have the 



following partition function for Cooper pairs: 

Z= fvpiVp2V9iV92VAe-'^, (Bl) 
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where p„ and On are the density and phase fluctuation 
of the n— th layer Cooper pair field, respectively, A and 
Aext are the fluctuating and external part of the elec- 
tromagnetic field, respectively. The intralayer Coulomb 
interaction Vi{r) = (2e)^/r (whose 2d Fourier transform 
would be 27r(2e)^/g), and the interlayer Coulomb interac- 
tion Ve{r) = {2e)'^ / y/r^ + a^ (whose 2d Fourier transform 
is 2TT{2e)'^ /qe^'^"'). ps is the superfluid phase stiffness of 
each layer. 

Similar to the single layer case in Appendix A, we can 
again introduce Hubbard-Stratonavich fields ji,2, split 
0's into smooth parts 9s and vortex parts 9y, integrate 
out 9s and A, and obtain 
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The difference from the single layer case is that now the 
continuity constraint is solved by introducing two new 
gauge fields a^^ = (ao,a) and /3^ = (/3o,/3) such that 
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Denoting the electric field and the magnetic field associ- 
ated with a^{P^) are Ci and &i (e2 and 62)7 respectively, 
we have 

ei = miji +.12) X z, 61 = rjic^i{pi + P2) 

62 = miJl - h) X ^: &2 = ??2C*2(/5l - ^2)- (B4) 

To fix 771,2 and the "speeds of light" c,i.2, we require 

47r ^Ps \ q J Ps q 

— {bj + bl) = -piV^pi + -p2ViP2 + PlVeP2, 

thus for n = 1,2, 
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g + gc(l-(-l)"e-9'^)' 



(B6) 



Using (B4) and (A15), we can again integrate by parts 
and express the partition function Z as 



Z = / VaVI3Ve^{De^2e 



~s 



(B7) 



where 



u),q 



S == ^i<-ielao + e2l3o) ipvi 



Bexi 
$0 



z(e*Q;o - 62^0) Pv2 - 



Be 



$n 



+ «Jd1 • (ei he2 )+«Jt,2-(ei e^ ) 

C*l C»2 C,i C,t2 



and for n = 1,2, the dual "charges" of the vortices are 



* ^ Trfi ^ / q 



(B9) 



When a (number) current bias j'l is applied in layer 1, 
the force on a vortex in this layer is 

F ^elx |ei| +ej X |e2| = 6^771 |Ji| +e^?72|Ji| 

and the force on a vortex in the other layer is 

F = ei X |ei| -- 62 X [eal = e*im\ji\ - e2V2\ji\ 
= 0, 



as expected. 

Again, introducing vortex fields V'^i and tpy2 for each 
layer and making the action explicitly gauge-invariant, 
we can write the action as in 



^ = E E 



q,LU 



2to„ 



+ Spyn {-ihuj(f)n - ie^ao + (-l)"2e2^o)] 



,1/2 2 2n / " 



4^"o + 4^^" 



+ ^(-^-^^2^ 4 



47r 



(BIO) 



Integrating out ag and /?o , one obtains the intralayer vor- 
tex interaction potential 



U.{q) 



*o9c 



q + qc 



and interlayer vortex interaction potential 



Ue{q) = - 



q + qc 



'u,.. 



(Bll) 



(B12) 



Which concludes the field-theory derivation of the inter- 
action potential. 



Appendix C: Classical derivation of the vortex 
interaction potential 

In this appendix, we present an alternative way of de- 
riving the vortex interaction potential between two vor- 
tices in a single superconducting thin film and in bilayer 
thin films. 

First, consider the current and electromagnetic field 
configuration of a single vortex at r = in a single su- 
perconducting thin film with thickness d located at z = 0. 
Combining the expression for the 3d current density of 
the vortex 



J = 



(^0- 



47rA2 \2'Kr 



e- A] 6{z)d 



(CI) 



where d is the thickness, and the Maxwell's equation, we 
have 



Next, we Fourier transform both sides of Eqn. (C2) 



(C2) 



^m k.) ^ ^-^^ (%, . = 0) - 1.;) , (C3) 

where q is the 2d wave vector, kz is the wave vec- 
tor in z— direction, and 6q is the azimuthal unit vector 
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in (7— space. Defining the inverse 2d screening length 
9c — d/{2\^) and integrating both sides j_ dk^, one 
obtains 



A{q,z^O) = — —0q 

q + qc iq 



From (CI), we have 



JW = — TT-^q- 

q + qc Sttz 



(C4) 



(C5) 



Now, we calculate the interaction potential between two 
vortices in a single superconducting thin film. The first 
vortex is located at r = 0, whose current distribution is 
given by (C5): 



Uq) 



qc C$0 , 



q + qc 27r« 
The second one is located at R away from the origin 

Mo) = J d\Ur)e-^^~'' ^ J d\Ur + R)e-'^^^ 

Their interaction potential is given by 
2n f d^q / 1 1 



(C6) 



(C7) 



UiR) = -J 



- + -]jii-q)j2{q}, (C8) 



c2 J (27r)2 yq 

where the first term is the kinetic energy contribution, 
while the second the term is from the magnetic energy 
B^ term. Using (C6) and (C7), we have 

c^ J {2t:Y \qc qj 



(2^)2 27r qiq + qc) 

where the vortex interaction potential 

$o9c 1 



(C9) 



U{q) 



27r q(q + qc 



(CIO) 



is exactly the same as what we obtained earlier in Ap- 
pendix A with field theory formalism. 

For the case of bilayer thin films with interlayer sepa- 
ration a, we can proceed in the same way. But there is 
one subtlety in that case. A vortex in layer 1, charac- 
terized by a phase singularity in layer 1, will also induce 
a circulating screening current in layer 2. Suppose the 
two identical layers are located at z = and z — —a, 
respectively, the one-vortex configuration is given by 



Ji 



Ji 



(^ 



47rA2 \27rr 
c 



e-A{z = Q)] 5{z)d, 



47rA2 



V'A- 



Att 



A{z = —a) ) 5{z + a)d. 



31+Ji 



(Cll) 




FIG. 5: The setup for calculating vortex interlayer interaction 
potential Us- A phase singularity in layer 1 leads to current 
ji and j'l in laye 1 and 2, respectively, and similarly a phase 
singularity in layer 2 leads to current J2 and J2 in layer 2 and 
1, respectively. 



Performing Fourier transform, one obtains 
2qc 



A{q,k,) = 



$0 



X I —e„ - A(q, z = 0) - e'^^^Aiq, z = -a 
iq 

Integrating over kz, one obtains two equations for 
A{q,z — 0) and A{q,z — —a), whose solution is given 

by 



A{q, z = 0) 



gc[g + gc(l-e-2'?")] 
[q + qc? 



q2(,-2qa 



£0 

iq 



-Ba 



Kq-, z^ -a) ^ 



Qcqe 



X —9„. 



(C12) 



[q + qc) ^ gcC ■"*" iq 

Thus, one can obtain ji and j( from (Cll) 

qdq + qc) c^o, 



ji 



(q + qc) 



fi = 



{q + qc 



2_q2^-2qa 27rj^''' 
2 _ ^2g-29a 27ri "?■ 



(C13) 



Next, one put in the currents j2 and J2 of another 
vortex either in the same layer or the other layer, and 
calculate the intralayer and interlayer vortex interaction 
potential Ui and Ue in the same way as we did for the 
single layer case. For example, to calculate the vortex 
interlayer interaction Ug, we put in another vortex with 
its core at the second layer, and it has a current J2 in 
the second layer, and a circulating screening current j'2 
in the first layer (see FIG. 5). Thus, 

(J1J2+J2J1) 

(C14) 



Ue{R) 



27r 

'^ J (27r) 

g-qa 



-{jih 



-j;j2) 



The final results are exactly the same as what we found 
in the field theory formalism in Sec. II A and Appendix 



B: 



U^{q) 
Uc{q) 



^O^Zc 



27r 9(g2 + 2(7c<? + g2(l-e-29«))' 



'U,.. 



(C15) 
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Appendix D: Classical hard-disc liquid description of 
the vortex metal phase 

As explained in Sec. II B, we expect that our results 
for the vortex drag do not depend sensitively on the mi- 
croscopic model we use for the vortices. In Sec. II B we 
used the fermionic vortex response function to determine 
the drag resistance in the intermediate metallic regime. 
Here we demonstrate the robustness of this result by re- 
producing the drag resistance results while modeling the 
vortex liquid in this regime as a classical hard-disc liquid. 

The density response function x(^i z) fo^ ^ liquid of 
hard-core disks in the hydrodynamical limit is^^'S^''^^ 



Xik,z) = x(fc) +iij;C{k,z), 



(Dl) 



where z is the frequency, T is the temperature, x{k) is 
the static compressibility, and 



C(fc,z) ^iTx{k) 



lz + iP(r + D(7-l)) 

7 



1 
1 - - 

7 



z^ — c^fc^ 
1 



zk'^T 



(D2) 



z + ik'^D _ 

showing a diffusive mode with weight 1 — - , and a prop- 
agating mode with velocity c, weight I/7 and life time 
l/(rfc2). Thus 



X{k,z) 



1 



Dk^ 



1 c^k^ - 
7 c^k"^ 



izDk^{j - 1) 



iTk^z 
(D3) 



X{k) V iJ Dk^ 
which satisfies the defining property of x'- 

X{k) ^\im x{k,z). (D4) 

1 is the constant volume specific 



Here, 7 = Cp/Cy, C^, 
heat, and 



Cp^Cy+TxTpv/n 



(D5) 



is the constant pressure specific heat, where n is the vor- 
tex density, xt = -^ linifc^o S(k) is the isothermal com- 
pressibility, and S{k) is the structure factor of the vortex 
liquid; j3v = n{l + y), where y = ^na'^g{(j), 



9{<y) 



1 - 7C/16 CV64 



{i-C? [i-CY 



(D6) 



c 



is the packing fraction, and a is the diameter 



of the hard-disc vortex which we take to be the core 
size of the vortex, which in turn is approximately 
superconducting coherence length ^ lOnm. 

In addition, T — a[ ^^^ ) + b, and the diffusion coeffi- 



cient D — -, where 



7' 



a^^ + ^l + iy/Afvl 



4 z/ 

,2 /Q I „,2/ 



(D7) 



b = 3i^cr'/8 + v^{l + y/2y/iy, 



V = 2y/TTnag{a)vQ is called the Enskog collision fre- 
quency, and the thermal velocity wq = y „; ?7i is the 
vortex mass. Finally, the speed of sound is 






(D8) 



The static compressibility x(fe) is related to the struc- 
ture factor S{k) (strictly speakly, the Ursell function'^^) 

by 



x(fc) = ^5(fc), 



(D9) 



and the structure factor S{k) of a hard disk liquid is 
determined by following the so-called Percus-Yevick ap- 
proximation of Ref 72,73: 



where 



S{k) = l/{l-nh{k)), 



h{k) = 27r / dRRJa{kR)h{R), 



^ ^ ' " R> 1 



0, 



(DIO) 



(Dll) 



(D12) 



Here, ^^ ~ 7r/16, C 



is the packing fraction. 



Hl)^^^'-'^^'-J^--f-^'-'^\ (D13) 



2{a- 13) 



h{0) = h{l) - I3h{l)^ , 



CSjR = 1) 

2/XD 



a = 2CA, 



(D14) 



(D15) 



A^ — (-\ f' dzz\l-z^)^^^,h = l + C, (D16) 
MD \aj Jo 



.X 1 I faR 

S[R] = - < arcsin —— 
a V 2 



d.R 



1- 



5.R 



-,1/2^ 



(D17) 



Putting these formulae together, we can compute the 
vortex density response function in (Dl) and insert it into 
the drag resistance formula (15). The drag resistance 
is shown in FIG. 6. One can see that it is remarkably 
close to our results obtained in Sec. (HB), and thereby 
demonstrating that the scale of the drag resistance in the 
metallic regime is mainly set by the factors dR/dB and 
is not sensitive to the statistics of the vortex particles. 
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B(T) 



FIG. 6: Drag resistance in the vortex paradigm at T = 0.07K, 
with the metalhc phase modeled as classical hard-disc liquid. 
Everything else is the same those in FIG. 2. 



Appendix E: Coulomb Drag for disordered electron 

glass 

In this section, we calculate the drag resistance due 
to Coulomb interaction between two disordered electron 
glasses with finite thickness. This calculation is related 
to the work of Ref. 74, but in our case the screening 
of the interlayer Coulomb interaction is important (see 
below), and we take into account the effect of finite film 
thickness. 

The general formula for Coulomb drag resistance in d 
dimensions is^"^'^'* 



11 
Pd 



1 1 



e2 27rn2T Vl ^^ 

k 



^'^^ / — rWl^l^IniXiImX2- 



sinh 



2T 



(El) 

For the quasi-2d film we are considering, we can break the 
wavevector summation into two summations: one over 
fc^, another over the 2d wavevector q. The k^ summation 
is dominated by the term with k^ — Q component, which 
physically corresponds to the configuration with constant 
density along z-direction. In this case, we can use the 
quasi-2d form of the intralyer and interlayer Coulomb 
interaction potentials 



C/,(q,fc, -0) 



2TTe^d 



Ue{q,k,=0) 



2TTe^d 



-qa 



where d is the film thickness, and a is the center-to- 
center layer separation. The real and imaginary parts 
of the density response function for a localized electron 
gas is 



75-77 



(g2 + kDu^S, 



2^f2 

2 I l,2\, ,ci 



2/-2 



Imx(<f, fcz = 0,a;) = v- 



D 



.2,.,(r4 



q'ojC 

= V 

fc^=o D 



where v is the 3d density of states at the Fermi energy, 
and ^ is the localization length, and D is the diffusion 
constant in the conducting phase. The above expression 
is valid so long as Imx ^ Rex, which is straightforward 
to verify in our case recalling that to is cut off by the 
temperature T in (El). 

Thus, in the screened interlayer interaction we can ne- 
glect Imx compared to Rex: 



U ■ 



Ke-«" 



(1 + C^»Xi)(l + ^^X2) - (^7^e-««Xl)(C^^e-9«X2) 

^ (E2) 



2Ui Re xi Re X2 sinh((7a) ' 

where in the last line we have made an approximation 
that t/jRex >• 1, i-e., 



ga > xo = 



v^^2TTe^d 



(E3) 



We have verified that the contribution from < go < xq 
is negligible compared to that from qa > Xq. Therefore, 



-Kd — — — 
d 



87r2(nd)2Te2 



duj 



q^dq 



Imxilmx2 



sinh^ — •^'"^2('T3 — ^2/o — \2„:„u2/ 



2T 



4C/2(Rexi)2(Rex2)2sinh"(ga) 
T2 



l2^'K^he^(ndaY{Die'^dv){D2e^dv) 
xdx f°° x^dx 



xa sinh X Jo sinh (x/2) 

r2 

128'K^he^{nda)^{Die^diy){D2e^diy) '"'' 2xo 3 



1 4^2 



T2 



IoSt;— 



96Tr^he^{nda)^(Die^diy){D2e^diy) ^ 2xo ' 

Note that 

1 



Deadly ^ 



we have 



i?' 



96TT^he^{nda)^ 2xo 

= ^^l^f^-Vlog^. 
967r2 h/e^ ye^nda J 2xq 



(E4) 



(E5) 



Since D is the diffusion constant in the conducting phase, 
R in the above expression should also be the resistance 
of the conducting phase. Thus this expression gives a 
slight overestimate of the drag resistance in the percola- 
tion paradigm if we use the value of Rnn of the insulating 
phase for simplicity. 

Note that our derivation relied on momentum sum- 
mations. There are concerns that such an approach, al- 
though quite common in the literature, is incorrect when 
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attempting to describe drag in strongly disordered sys- 
tems. For our purposes, the derivation based on Eq. El 
is sufficient; this issue is taken up, however, in Ref. 78. 



Appendix F: The absence of measurable drag effect 

associated with a genuine superconductor in a 

resistor networli 

In this section, we show that a genuine superconduct- 
ing hnk (i.e., without mobile vortices) has no measurable 
drag effect in a resistor network. 

K 



i?, 



Active Layer 



and aD,NS = v/Ri < co. From Kubo formula for the 
drag conductance, we expect that aD,SN = o'd,ns < cOj 
and hence for the case of i?i =0 and i?2 > we have 
V = o-d,snRi = 0. 

In contrast, the second drag effect is the drag current 
I2 in the presence of Rq , in which case he drag current at 
Rq — may or may not survive. In a large-size resistor 
network we are considering for the percolation picture, 
when we focus on the drag effect of one specific link R2, 
we can simplify the circuit of the passive layer to be of 
the form in FIG. 7, in which case Rq representing the 
rest of the circuit is almost always larger than 0. If the 
drag effect survives the presence of the nonzero Rq, it 
will manifest itself as the appearance of a non-zero drag 
emi Ed on i?2. To see this, first consider the case R2 > 0, 
and i?i can be either or > 0. I2 receives contribution 
from both Ohm's law and the drag effect: 



R2 



Passive Layer 



V, 



R. 



FIG. 7: The typical setup for a drag effect experiment: in the 
active layer, a driving current /i flows through a resistor Ri 
(normal or superconducting) with a voltage drop Vi — IiRi- 
In the passive layer, certain interaction effect takes place in a 
resistor R2 (normal or superconducting) , which may result in 
a drag current I2 and a voltage drop V2 across R2 . R2 is also 
connected to another resistor Rq, which can be of any value. 

FIG. 7 illustrates the typical setup for a drag effect 
experiment: in the active layer, a driving current /i flows 
through a resistor i?i (normal or superconducting) with 
a voltage drop Vi = IiRi. In the passive layer, certain 
interaction effects take place in a resistor R2 (normal or 
superconducting), which may result in a drag current I2 
and a voltage drop V2 across R2 . R2 is also connected to 
another resistor Rq, which might represent a voltmeter, 
an open circuit (_Ro = 00), or something else. 

When one talks about the drag effect, there are two 
different concepts one needs to distinguish. The first 
one is the "intrinsic" effect, which manifests itself by the 
appearance of a drag current Id in the passive layer if 
Ro — 0- Gcnerically, we have 



'2|_R,o=0 



vh 



(Fl) 



For example, for the case of i?i,-R2 > 0, i.e., both i?i 
and i?2 arc non-superconducting, /2|flo=o = <7dVi = 
(JdIiRi (e.g., Coulomg drag between two 2DEGs), thus 
rj = ubRi', for Ri = R2 — Q (superconductor), we have 
the Cooper pair version of the supercurrent drag effect 
Eqn. (32), thus 77 is finite in this case as well. For the 
case of i?i > (normal) and -R2 = (superconducting) , it 
would be unphysical to have rj — 00, thus we have 77 < 00 



R2 



Vh 



I2R0 



77/1 



thus 



/2 = 



{vR2)h 



Roh 



£d 



Rq + R2 -Ro + R2 Rq + R2 



(F2) 



(F3) 



where Sd — RdIi is the drag emf, and Ro — tjR2 is 
the drag resistance. If i?i =0 (superconducting) and 
i?2 > (normal), we argued earlier that 77 = 0, and thus 
Ed = Rd = and there is no drag effect. 

If i?2 = (superconductor), no matter if _Ri =0 (su- 
perconducting) or > (normal), it is straightforward to 
see from Kirchoff's Law that we have only one steady- 
state solution /2|_Ro>o ~ 0. More insight into this case 
can be gained by considering what happens in real time. 
Suppose at time t — 0, the drag effect takes place, a drag 
supercurrent l2{Ra = 0) starts to flow in the circuit. But 
due to the presence of the normal resistor _Ro, a voltage 
/2-R0 now exist on the supercondutor, which will crank 
up the phase winding of the superconductor and degrade 
the drag supercurrent, until a steady state is reached 
where the total supercurrent is zero. Thus, we see that 
for the case _R2 = and i?o > 0, there is no observable 
drag effect, i.e., /2Uo>o = 0, ££) = /2(i?2 + Ro) = 0, 
Rd = £d/Ii = 0, although there is nonzero "intrinsic" 
drag effect 77. 

We can also understand this result Rjj = for _R2 = 
by examining the expression Ru ~ 77i?2- For both the 
case of i?i = i?2 = and the case of i?i > and R2 = 0, 
we found earlier that 77 < 00, and thus the drag resistance 
Rd = riR2 a-nd the drag emf sd are for R2 = 0. 

In conclusion, we have shown that when connected 
with a nonzero resistor, as typically true in a resistor 
network, a genuine superconducting link has no measur- 
able drag effect at all, no matter whether it is vertically 
aligned with a normal link or another superconducting 
link. 
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